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^S) ■ We discuss a number of lattice fermion actions solving the Ginsparg- Wilson 

relation. We also consider short ranged approximate solutions. In partic- 
ular, we are interested in reducing the lattice artifacts, while avoiding (or 
suppressing) additive mass renormalization. In this context, we also arrive 
£^ , at a formulation of improved domain wall fermions. 

00 

Q\ ' I. The remnant chiral symmetry of fermionic fixed point actions 

It is a notorious problem to find a lattice regularization for fermions, 
q_|' which preserves chiral symmetry and avoids species doubling. The cele- 

brated Nielsen Ninomiya No-Go theorem excludes the reproduction of the 
chiral symmetry manifestly in a local lattice action of a single fermion, if 
some plausible assumptions are respected Q (here, locality means that the 
^ ■ inverse propagator is analytic in momentum space). However, the manifest 

chiral invariance of the lattice action is not badly needed to its full extent: 
with respect to important issues, it is sufficient to preserve just a remnant 
chiral symmetry, circumventing the Nielsen Ninomiya theorem. 

One way to tackle the problem uses the technique of block variable renor- 
malization group transformations (RGTs). Under such transformations, the 
partition functions and all expectation values - hence the physical contents 
of a theory - remain invariant Q. Usually one maps the theory from some 
fine lattice to a coarse one, increasing the lattice spacing by a factor n. An 
infinite number of iterations at infinite correlation length leads - for suitable 
RGT parameters - to a fixed point action (FPA) of the considered RGT. 
If we send the blocking factor n —* oo, and we refer to coarse lattice units 
from the beginning, then we map a continuum theory on the lattice, without 
altering its physical contents. This mapping, which generates the FPA in 
just one step, involves a functional integral over the continuum fields. 

For free fermions, such a blocking from the continuum can be performed 
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by the RGT 



e -S[%9] = / D ^ D ^ e -s[^] x (1) 

ex P {- J2 [** - / drU ^ x - r )^( r )} <*xb [% - [ dr'Tl(y - r')^(r')] }, 

where ^, ^ are the continuum fields, ^, the lattice fields, s is the con- 
tinuum action, and S the fixed point lattice action. The matrix a and 
the convolution function II specify the RGT (the latter is peaked around 
and normalized as J drIL(r) = 1). For a massless fermion, we obtain in 
momentum space the fixed point propagator P, f|] 



, U 2 (p + 2irl) , , 



7T,7Tl d . 



It turns out that the FPA characterized by this propagator is free of fermion 
doubling [§, §. 

In the limit a = 0, the blocking of the fermion fields is implemented by 
5 functions. In this case, and generally for any RGT with {0,75} = 0, the 
lattice action is chirally invariant ({,} denotes the anticommutator). There 
is no contradiction with the Nielsen Ninomiya theorem, because in all these 
cases the action is non-local ||: G~ 1 (p) has poles at the corners of the 
Brillouin zone B. This is a non-locality of the type proposed by Rebbi ||. 
His fermion proposal was dismissed, however, because it implied zero axial 
anomaly |0j, and - as an extension of the Nielsen Ninomiya theorem - it 
has been shown that this is the case for a whole class of non-local fermions 
with poles in G~ 1 (p) However, a consistent blocking of the fermionic 
fields, of the free gauge fields, of the interaction term to the first order in 
the gauge coupling, and of the axial charge and current, does reproduce the 
axial anomaly correctly on the lattice ||. [] 

Now we turn to the case where the action (and the axial current) becomes 
local, so we now deal with {a, 75} 7^ 0. In this case, which is more interesting 
in view of practical applications, the RGT term is not chirally invariant, and 
therefore the chiral symmetry in its naive form is not manifest in the fixed 
point action. Hence the No-Go theorem does not apply, but the chiral 
symmetry is still represented correctly in the observables, due to the very 
nature of the RGT || . This is an elegant way to by-pass the No-Go theorem. 

Since integrating out the continuum fields requires a functional integral, 
the FPA - and more generally, the perfect action - can be made explicit in 
perturbation theory. Beyond that, also non-perturbative properties, which 
are known in the continuum, can be reproduced correctly on the lattice, but 

1 The reason why the generalized No-Go theorem does not apply is that the consistent 
perfect axial lattice current is also non-local, whereas it has been assumed to be local in 
the proof of Ref. fl. 
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the expressions for the lattice action seem to become somewhat symbolic. 
However, it is very interesting that the Atyiah-Singer index theorem still 
holds even for the classically perfect gauge action |H| (used together with 
a classically perfect topological charge fll|| ). Still, the action described in 
this way (in terms of classical inverse blocking by minimization) is somehow 
implicit, but certainly not symbolic any more. In the Schwinger model, it 
was indeed possible to confirm numerically the index theorem for a classically 



perfect action (to a good accuracy) [12] 



The fixed point propagator (y) obeys 

{G, 75 } = {«,75}. (3) 

In the case of a local - but not chirally symmetric - term a, this can be 
viewed as a remnant chiral symmetry of the lattice action. This relation has 
been postulated by Ginsparg and Wilson as the 'softest way' to break chiral 
invariance in the lattice action In this form, the remnant chiral symmetry 
is particularly transparent: it is inherent to the propagator up to a local 
violation. ^ In particular, it was demonstrated in Ref. || that the Ginsparg- 
Wilson relation (GWR), eq. (||) with local anticommutators, is sufficient to 
obtain the correct triangle anomaly, and also the soft pion theorems were 
expected to be reproduced correctly. Furthermore, this property is sufficient 



to obtain the correct anomaly of Tr{^G ) EQ, 14]. In addition, it has 



been shown that eq. (|3|) implies a continuous symmetry, ^ since the action 
is invariant under the substitution [14] 



^ -> (1 + e[l - Grafts) , * -» (1 + e 75 [l - aG' 1 }^ (4) 

to 0(e), as we see if we write the GWR in the form 

{ 75 , G- 1 } = G- 1 {75, ajG- 1 . (5) 

All this sheds light on the "miraculous way" the FPA finds to circumvent 
the No-Go theorem. 

Actually, the transformation (|j) is a generalization of the case a = 1/2 



considered in Ref. [14], to any a permissible in the GWR. However, in 
the following we also put special emphasis on a = 1/2. If we perform the 
blocking (jl]) with the standard block average scheme (n(u) = 1 if \u^\ < 
1/2, n = 1 . . . d, and n(u) = otherwise), then the choice a = 1/2 (or 
a = —1/2) does in fact optimize the locality of the free FPA, which we 

In contrast, the usual chiral symmetry breaking terms like mass or a Wilson term 
correspond to a local violation in the inverse propagator. There the symmetry breaking 
in the propagator becomes non-local. The same happens if we block a massive fermion 
from the continuum (the mass is added in the denominator of the propagator (Q) p3|). 

3 Blocking a SUSY theory from the continuum could reveal an analogous remnant lattice 
supersymmetry (but the resulting action differs from the one proposed in Ref. [^5[ ). 
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write in coordinate space as 



S[<M]= *Mrh„ + \(r)]V x+r . (6) 

x,r&% d 

For this choice, the couplings p^(r), A(r) are restricted to nearest neighbors 
in d = 1, and in d > 1 their exponential decay is extremely fast. ^ The 
numeric values of the leading couplings in d = 4 are given in Ref . |l3| . 

If we consider the form (^) as a general ansatz for the lattice action, 
normalized so that J2r r ^P^( r ) = ( n o sum over /i) and J2r ^( r ) = 
(which means that the fermion is massless), then we see that it is hardly 
possible to find an ultralocal solution of the GWR (by "ultralocal" we mean 
that the couplings are restricted to a finite range). In general, A may have 
a Dirac structure, but in the following we assume its form to be scalar (in 
Dirac space). Thus the GWR (||) with a = const, simplifies to 

A(r) = a [Kv)H r -y)- p»(y)p»( r - v)\ ■ (7) 

ye% d 

For an ultralocal action, the right-hand side extends at least to twice the 
range of the left-hand side (it extends even further if a is momentum de- 
pendent) and it is very unlikely that the latter can achieve cancelations 
everywhere beyond the range of A, match with the left-hand side at each 
point inside that range, and obey the correct normalization. [] All this tends 
to strongly over-determine the degrees of freedom in an ultralocal action. 
This suggests that in the most local exact 2d and 4d solutions, the couplings 
in G _1 (r) decay exponentially in \r\. [] 

Considering eq. (^), one is tempted invent new solutions by hand, such 

as 

G{p) = — + a , G (p)= f f 2 _ 2 8) 

where p^ is some lattice version of p„, and a is a scalar. But for instance 
the obvious example p^ = sinp^, which corresponds in fact to a local action 
for any real a, suffers from fermion doubling. Now one has to remove the 
doublers, maintaining the locality of {^,75} (and preferably also of G _1 ), 
and it is not trivial to find sensible alternatives to the fixed point propagator 



4 For instance, on the 4-space diagonal the decay behavior is p^(n, n, n, n) oc exp(— 4.94- 
n), A(n, n, n, n) oc exp(— 4.97 • n). 

5 Summation over r leads to the requirement Pn(y) = 0. This is guaranteed for the 
sensible assumption that p M is odd in y^, which holds for the FPA in the block average 
scheme with an even Dirac scalar a. There, p M is in addition even in all ?/„, v ^ fj,, and A 
is even in all directions. 

6 The simplification a la eq. ([?]) for a = ±1/2 is solved in d = 1 by the Wilson action 
with Wilson parameter rw = ±1. 
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r 


75 {G- 1 (r), 75 } 


Is^G-^bsG-^r-x) 


d = 4 






(0000) 


3.70544 


3.70544 


(0001) 


-0.12152 


-0.11735 


(0011) 


-0.06007 


-0.06014 


(0111) 


-0.03194 


-0.03356 


(1111) 


-0.01685 


-0.01805 


(0002) 





-0.00417 


(0012) 





0.00020 


(0112) 





0.00069 


(1112) 





0.00047 


(0022) 





-0.00034 


(0122) 





-0.00013 


(1122) 





-0.00008 


(0222) 





-0.00006 


(1222) 





-0.00005 


(2222) 





-0.00002 


d= 2 






(00) 


2.97909 


2.97909 


(01) 


-0.48955 


-0.48631 


(11) 


-0.25523 


-0.26035 


(02) 





-0.00324 


(12) 





-0.00350 


(22) 





-0.00188 



Table 1: The two sides of the Ginsparg- Wilson relation (]5|) with a = 1/2 for 
the truncated block average FPA ("hypercube fermion") in d = 4 and in d = 
2. For all vectors r, which do not occur in the table, both quantities vanish. 
The exact agreement at zero distance is a consequence of the truncation by 
periodic boundary conditions. 



However, if we can not achieve ultralocality, in practice a truncation of 
the couplings to a short range is needed, which causes a violation of the 
GWR. Therefore, what one should aim at is a truncation that keeps this 
violation small, and here the optimization of locality helps. The trunca- 
tion can be carried out in an elegant way by means of periodic boundary 
conditions, which keep the normalization exact. This truncation to cou- 
plings in a unit hypercube - for the block average FPA with a = 1/2 - 
has been presented in Ref. 17] ("hypercube fermion"), and its spectral and 
thermodynamic properties are drastically improved compared to the Wilson 
fermion [17, Gil]. Table [l] shows that also the GWR is approximated very 
well by this truncated perfect fermion. 

This provides hope that the desired properties related to the remnant 
chiral symmetry are realized to a good approximation for the "hypercube 
fermion" (HF), in particular the absence of additive mass renormalization. 
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That property is exact for the perfect and for the classically perfect action 



[19|. While the index theorem can be considered as classical, it is especially 
remarkable that also a quantum property like the stability of the chiral 
limit under renormalization is fulfilled by the classically perfect action. In- 
tuitively, it appears that in the chiral limit the classically perfect actions 
also displays properties of quantum perfection, and it is plausible that the 
continuous remnant chiral symmetry (||) protects the zero bare mass from 
renormalization, similar to the remnant chiral symmetry U(l) ® U(l) of the 
staggered fermions. Here we have the additional virtue that the number of 
flavors is arbitrary (since doubling is avoided). 

As a further virtue of the classically perfect action in the chiral limit, 
there are no exceptional configurations as shown in Ref. (in the sense 
they are understood there). However, this paper refers to a fixed point ac- 
tion, which is defined implicitly by minimization. For practical applications, 
we need a very simple gauging prescription (involving only a few short lattice 
paths), and the requirement to preserve the GWR to a good approximation 
could guide the construction of such a good but simple gauging. We com- 
ment on this concept in the Appendix. 

II. A new class of solutions of the GWR 

A different solution of the GWR arose from the so-called overlap formal- 



ism [ 20 [ . To turn it the other way round, for any solution of the GWR, the 
fermion determinant factorizes in a way compatible with a vacuum overlap 
plf . Also this consideration focused on a = 1/2, and it is interesting that 
an entirely different approach singles out the same a as the optimally local 
block average fixed point fermion. 

We re-derive the solution of Ref. p0| ] and reveal its minimal assumption, 
which allows for a broad generalization. First we define 

V = 1 - G- 1 , (9) 

and eq. (||) with a = 1/2 is equivalent to 

75^75 = V- 1 . (10) 

For the case of one flavor, the solution of Ref. |2(|] reads 

V = X{X ] X)- 1 ' 2 , X = 1 - D w , (11) 

where Dw is the standard Wilson-Dirac operator. 

Depending on the gauge background, there is a danger of a zero eigen- 
value of X, so that the expression for V is undefined, but following Ref. [^] 
we assume that this can be ignored statistically. Note that whenever V is 
well-defined, then it is unitary. Hence all we need for eq. (pi]) to hold, is 
the property 

7 5J D W 75 = D w , (12) 
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since this behavior is inherited by X and V (and by G~ l ). This leaves a lot 
of freedom for generalizations. The minimal condition is 75X75 = X\ but 
we stay with the form 

X = 1 - D (13) 

and note that one may insert many lattice Dirac operators for D. If we 
make also here the ansatz D{r) = /0«(r)7„ + A(r) - with the normalization 
mentioned before, and with a scalar form of A - then this property always 
holds. 

Now the question arises, what generalization is useful. Locality is a 
criterion, but even if D is ultralocal - as in the case D = D\y - then the 
decay of the couplings in G _1 (r) is still exponential (in d > 1). For Dyy this 
can be seen from 

G-\p) = 1 - [l - isinjyy,, - ^p 2 ] x (14) 



2 



1 ^ rl 



1 + (i - r w )f +ye^ + -fHpipi 

(J. \i>V 



-1/2 



which simplifies for the standard Wilson parameter rw = 1 and for ry/ = — 1 
(the value that the domain wall literature deals with). Locality is realized 
for any choice of rw [] 

Returning to our general ansatz for D, we obtain 

G-\p) = 1 - [1 - ppip)^ - A(p)][l - 2A(p) - p 2 ( P ) + X 2 (p)]- 1/2 , (15) 

and together with eq. (0) this leads to an amusing observation: if we insert 
an operator D, which obeys the GWR itself, then X becomes unitary, and 
we end up with 

G- 1 = D . (16) 

For example, the inverse fixed point propagator reproduces itself identically, 
if we insert it for D. If D violates the GWR, then it gets "GWR corrected". 
This correction cannot cure all diseases, however: if D is non-local or plagued 
by doubling, then the same is true for G _1 . On the other hand, also massive 
fermions can be GWR corrected. 

As a further criterion, one should aim at small lattice artifacts. This 
suggests to insert a choice for D, which has good properties in view of 
approximate rotation invariance etc., because such properties are essentially 
inherited by G~ l , as the expansion (|l~5|) shows. Hence, as an experiment we 
insert the HF for D. Its small GWR violations (see Table |l|) are corrected 
in but an exponential tail of couplings is added, so that we obtain 

something similar to the original fixed point propagator. Of course, this 
G^ 1 is not fully perfect, but let's assume that for the aspects of interest 



7 For rw = 1/d the form of G 1 (r) simplifies to <5 r ,o+ even-odd-couplings, but the latter 
decay rather slowly. I thank M. Peardon for this remark. 
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FPA 


HF 

GWR 


r w = 0.85 
GWR 


rw = 1 
-> GWR 


rw = 1.15 
-> GWR 


r w = -1 
-> GWR 


d = A 














r P 


1.635 


1.519 


2.688 


2.530 


2.537 


2.424 


r\ 


1.187 


1.109 


1.819 


1.708 


1.810 


2.784 


d = 2 














r P 


1.268 


1.198 


1.981 


1.816 


1.844 


1.946 


r\ 


0.871 


0.844 


1.286 


1.177 


1.206 


2.330 



Table 2: The locality of the vector and the scalar part of the FPA, 
of the GWR corrected version of the "hypercube fermion", and of the 
GWR corrected Wilson fermion with various Wilson parameters r\y. (The 
GWR corrected actions are not normalized, because otherwise they vio- 
late the GWR again.) As a characteristic quantity, we measure r 2 p = 
(Ex \Pn(x)\x 2 )/(Ex IPm( x )D' and analogously r x . 



only the GWR matters. Then the action constructed in this way looks 
useful, because it is more local than the original FPA (and much more local 
than the actions obtained from D = Dyy), see Table ^. Requiring only the 
GWR means to relax the condition of perfection a little, and we make use of 
that by further improving locality. For comparison, the leading couplings 
for some of the fermion actions discussed here are listed in Table ^. 

In order to generalize this class of GWR solutions to any (non-trivial) 
local Dirac scalar a, still assuming 75-D75 = D* , we write 

X = 1 - 2^ D y/a , G^ 1 = — ?= [l - X(X^X)- 1 / 2 } -L • (17) 

2v« L J 



Again, G^ 1 fulfills the GWR, and if also D does so, then we end up with 
G _1 = D, since X is unitary. The latter further implies that the spectrum 
of any solution G , at a specific momentum p G B, is situated on a circle 
in the complex plane, 

<G-\p)) C ^-(l + e*) , (f £ [0, 2ir[. (18) 

The entire spectrum of G^ 1 is located between the circles given by the 
minimum and the maximum of a(p). f\ For our preferred choice a = 1/2, 
the entire spectrum is located on one unit circle. In Figs. |] and [2] we 
illustrate that this property is still approximated very well for the FPA 
after truncation. We compare the HF and its GWR corrected version, which 
is truncated again in coordinate space, given in the second and in the last 

8 The same geometric structure was obtained in Ref. |l(| for the FPA with constant a, 
but iteration of overlapping blocking functions IT. 
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FPA 


HF 


HF 
-> GWR 


rw = 1 
-» GWR 


r w = -1 
-> GWR 


HF -> GWR 
trunc. again 


d= 4 














Pi (1000) 


0.1400 


0.13685 


0.1363 


0.1668 


0.1023 


0.13633 


pi(HOO) 


0.0314 


0.03208 


0.0318 


0.0239 


0.0105 


0.03185 


pi(1110) 


0.0096 


0.01106 


0.0110 


0.0050 


0.0021 


0.01108 


Pl(llll) 


0.0033 


0.00475 


0.0048 


0.0007 


0.0006 


0.00488 


pi (2000) 


0.0040 





0.0002 


0.0238 


0.0107 





Pi (1200) 


0.0015 





-0.0001 


0.0053 


0.0017 





A(OOOO) 


1.8530 


1.85272 


1.8534 


1.8334 


0.0471 


1.85313 


A(1000) 


-0.0616 


-0.06076 


-0.0626 


-0.0560 


0.0093 


-0.06289 


A(1100) 


-0.0290 


-0.03004 


-0.0300 


-0.0241 


0.0027 


-0.03024 


A(1110) 


-0.0146 


-0.01597 


-0.0152 


-0.0115 


0.0010 


-0.01545 


A(llll) 


-0.0077 


-0.00843 


-0.0079 


-0.0057 


0.0004 


-0.00813 


A(2000) 


0.0010 





0.0019 


-0.0002 


-0.0038 





A(1200) 


-0.0007 





-0.0002 


-0.0029 


-0.0005 





d= 2 














Pi(10) 


0.3154 


0.30939 


0.3089 


0.3330 


0.1680 


0.30907 


Pi(H) 


0.0870 


0.09531 


0.0953 


0.0580 


0.0273 


0.09540 


pi (20) 


0.0076 





0.0000 


0.0459 


0.0288 





Pi(12) 


0.0048 





0.0002 


0.0166 


0.00068 





A(00) 


1.4903 


1.48954 


1.4899 


1.4912 


0.0651 


1.48955 


A(10) 


-0.2441 


-0.24477 


-0.2464 


-0.2335 


0.0173 


-0.24683 


A(ll) 


-0.1246 


-0.12761 


-0.1252 


-0.1167 


0.0061 


-0.12556 


A(20) 


-0.0009 





0.0015 


-0.0128 


-0.0103 





A(12) 


-0.0032 





-0.0017 


-0.0108 


-0.0028 






Table 3: The leading couplings of some actions discussed in the text: the 
FPA, the "hypercube fermion", the "GWR corrected hypercube fermion" 
and the "GWR corrected Wilson fermion" with rw = 1 (standard) and 
rw = — 1 (standard for domain wall fermions). All these actions have the 
symmetries described in footnote 5. The similarity in the first three cases 
is not surprising, but it is amazing that also the fourth example looks quite 
similar. For the Wilson fermion we could also vary a according to eq. (|T7D- 
If we truncate the "GWR corrected hypercube fermion" again periodically, 
then we reproduce the "hypercube fermion" identically, hence we truncate 
in coordinate space this time and then correct the normalization by hand; 
this causes less alteration and yields the last column. (An iteration of GWR 
correction and truncation does hardly modify the couplings any further.) 
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column of Table |3[ Generally, the eigenvalues are given in d = 2 by £1,2 (p) = 
X(p) ± i\J p' 2 (p), and in d = 4 by the roots of 

£ 4 Cp) - 4X( P )e 3 (p) + [6A 2 (p) + 2p\p)}e 2 {p) 
-AX(p)[X\p) + p 2 (p)]e(p) + [X 2 (p) + ^(p)] 2 . 

While the progress due to GWR correction is evident in d = 2 (Fig. |l|), the 
difference is not easily visible in d = 4 (Fig. ^). However, if we measure the 
"mean deviation" from the unit circle by 

2 d/2 

52 = J^r L dp lw E ^ - x l 2 ] > ( 19 ) 

then we observe that 5 decreases from 1.09- 10~ 2 for the HF to 0.75- 10~ 2 for 
its GWR corrected and truncated modification. This shows that the latter 
solves the GWR to an even better approximation. 

For the free fermion, p = always yields £j = 0, i = 1 . . . 2 rf//2 , and the 
whole sector e ~ is hardly affected by the truncation (see Figs. 1 and 2). If 
we add a gauge interaction, then this sector is crucial for the index theorem 
and for the absence of "exceptional configurations" , hence its good quality is 
very important. Truncation effects are rather visible at the opposite sector 
e ~ 2, which corresponds to large momenta. 

In a direct application of HFs, it turned out that even using a strongly 
simplified gauging (instead of the consistently (classically) perfect incorpo- 
ration of the gauge field), for instance the meson dispersion is drastically 
improved Jl7| , |22"| . However, this application suffers from a strong mass 
renormalization, which is a problem with respect to practical issues. In par- 
ticular, for the minimal gauging by hand along shortest lattice paths, the 



"pion" mass is renormalized from 0.0 to 3.0 at f3 = 5 [17|. Thus one can 
hardly trust the couplings in general, and in addition tuning to the critical 
bare mass is required. Its sign is opposite to the sign of a, which means that 
we are led to a region, which is unfavorable for locality. Since the truncation 
effects in the free fermionic FPA are small, we expect that this renormal- 
ization could be almost avoided by using a quasi (classically) perfect gauge 
interaction, which is, however, very difficult to realize in d = 4. If we gauge 
by hand, trying to suppress the additive mass renormalization (as in Ref. 
p3[), then one should use the most local free fermion couplings, and here 
the "GWR corrected hypercube fermion" is in business. For the idea to use 
the GWR as a guide-line to construct a better gauging by hand, we refer 
again to the appendix. 

III. Improved domain wall fermions 

So far we have considered the construction of a potentially interesting 
GWR corrected propagator for the direct use. However, the way people 
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think about simulating overlap fermions uses directly D in the above lan- 
guage. In that context, the quality of the locality of G' 1 is not directly 
relevant, hence one could choose D = D\y and then gauging is not problem- 
atic either. On the other hand, one runs into some trouble with 'exceptional 
configurations' (now in the sense that X has a zero eigenvalue). f\ More im- 
portantly, an extra dimension has to be introduced, which can be interpreted 
as a fifth direction with L s sites, or simply as L s flavors. The meaning of 
the square root, that the fermion determinant is divided by, is to subtract 
the contributions of heavy modes, generated by the extra dimension. In 
practice one has to work with auxiliary boson fields [24]. 



For practical purposes, the domain wall fermion formulation by Y. Shamir 
seems most suitable [25], hence we want to discuss possible improvements 



in that framework. In terms of vector-like gauge theories, first simulations 



have been performed for the Schwinger model pq , 27] and for QCD [28]. 
Again, an extra dimension is required, and the limit L s — ► oo has attractive 
features: absence of additive mass renormalization (supported by a one loop 
calculation f2j|), for the measurement of matrix elements no fine tuning is 
needed either |p8| , and there are only 0(a 2 ) artifacts 0. In practice L s must 
be finite, so that these properties are not exact anymore, but one expects for 
instance the 0(a) artifacts to be suppressed exponentially with L s . Also the 
chiral limit appears to be stable to a decent approximation down to L s ~ 10 
[p8|. Questions of the L s truncation have been discussed on the theoretical 



level in Refs. [25, 3C[. We note that Shamir's formulation allows to work 



with L s /2 only, in contrast to Kaplan's original proposal [31]. 

The action, which has been used so far, has the following form |p2|| : in 
Euclidean space - at a fixed extra coordinate s - it is the usual Wilson action 
with gauge fields on the links, with ryy = — 1 and a mass M. In the extra 
dimension, it is the free massless Wilson action with ry/ )S = — 1 (the kinetic 
term with 75), and with one additional term 

D a s d s> = \ [(1 " 75)^,o^,l s -i + (1 + 76)^,0^.-1] (20) 

(we refer to one flavor). The bare quark mass is m q = M(M — 2)(1 + a), 
but due to renormalization, there is a tuning problem for M. On the other 
hand, one expects (and observes reasonably well p8| ) m\ oc (1 — a). 

Now we distinguish three possible strategies for the improvement of this 
type of domain wall fermions: 

1) The observation that certain artifacts of Dyy, for instance in the 
rotational invariance, essentially persist, motivates the use of the HF - or of 

9 The extent of this danger depends on rw'- for rw < 1 the smallest eigenvalue of 
X ] X in the free case is 1, but for larger rw it decreases down to a minimum of 0.288 at 
rw = 6.4, before it rises again to 1. 

10 This is related to the chiral symmetry at L s — > 00, which rules out all 0(a) operators. 
Again, there is a certain similarity with staggered fermions, but here we keep the whole 
SU(N) ® SU(N) symmetry, hence the correct number of Goldstone bosons is involved. 
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its GWR corrected and truncated modification - in the 4d (or 2d) Euclidean 
space. Then for example a simplified gauging would be less problematic, 
because the additive mass renormalization and other diseases are suppressed 
by L s . Still, such a simulation would be tedious, but so is the use of a 
properly gauged HF. If this works, then one cumulates all sort of advantages: 
very small lattice artifacts (doubly suppressed from the HF and from L s ), a 
small additive mass renormalization, an arbitrary number of flavors and the 
correct number of Goldstone bosons. It remains to be checked if the left- 
and the right-handed fermions can still be separated, see below. 

2) Instead one could improve in the 5-direction. If we assume L s peri- 
odicity, then its inverse propagator in the above Wilson-like form reads 

D W,5(P5) = i sinps75 + r -^-v\ + y ex P( iL sP575) . (21) 

If we block in the 5 direction only, the corresponding "FPA" is given by P 

where II5 is the blocking prescription in the 5-direction. If we choose again 
the block average, and we optimize 05 for locality, then we arrive at 

D^{p 5 ) = (l + a)(|) 2 [isinp 5 % ± (ip§ + «)] (23) 

u = e u - 1 , p 5 = -^j, j = . . . L s - 1. 



(1 + a)L s '° L s 

In a sense, there in not much to improve in one dimension, since the Wilson 
action (with rjy = ±1) is already perfect; we have just incorporated and 
L s . However, the s dependent defect, which is needed to separate the chiral 
fermions, has disappeared. A way to avoid this is the use of fixed boundary 



conditions in the s direction, similar to Ref. [33]. This study is in progress. 

3) Finally, one could improve directly the entire 5d action, allowing for 
diagonal couplings with respect to all directions. Since there are no gauge 
fields on the 5-links, this is perhaps not as disastrous as it first appears. 
With respect to gauging, there is no additional problem compared to the 
case 1), and we hope to suppress all discretization artifacts at once. Here, 
a conspiracy of M, u and L s affects the couplings in all directions. 

Also in this case, L s periodic boundary conditions destroy the defect in 
the s direction 0, hence a more sophisticated treatment is needed. 



FPA" needs inverted commas, because we are at finite L 3 , but let's assume that we 
start at a huge size n k L s , then perform k RGTs with block factor n, so we end up in L s . 
In the limit n k — > 00 we obtain the expression (0). 

12 This happens even if we split oT 1 into [a£ 1 (l — 7s)+a i j 1 (l+75)]/2, possibly combined 
with chirally specific functions ITl, II_r. 
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r 


pi(r) 


A(r) 




d = 4 


M = 1.7 


(0000) 





0.95860 


(1000) 


0.02606 


-0.01613 


(1100) 


0.00468 


-0.00476 


(1110) 


0.00123 


-0.00178 


(1111) 


0.00041 


-0.00076 




d= 2 


M = 0.9 


(00) 





1.14008 


(10) 


0.14008 


-0.10247 


(11) 


0.03278 


-0.04382 



Table 4: The truncated perfect "hypercube fermion" at the masses, which 
have been used in domain wall fermion simulations: M = 1.7 in d = 4, and 
M = 0.9 in d = 2. 



The defect in the s direction is needed to keep the chiral fermions apart, 
so we focus on variant 1), which leaves the couplings between the 4d layers 
untouched. Since in general a mass parameter M appears in the physical 
subspace, we insert a massive HF. This is constructed from the perfect action 
at mass M, again by truncation with periodic boundary conditions [17]. 
There are two choices for a which restrict the Id action to nearest neighbors, 
a = au = (e — M — 1)/M 2 > 0, and a = —a-u- In higher dimensions, 
the first (second) option optimizes the locality for positive (negative) M. 
□ In the following, we choose a = au- As a sound basis, the hypercubic 
couplings for typical masses, which have been used in simulations, M = 1.7 
in d = 4 (the approximate (quenched) critical value at (3 = 6 [^]) and 
M = 0.9 in d = 2 27|, are given in Table [|. Note, however, that the 
renormalization with an improved action could lead to a different critical 
value of M. 

Finally, we want to show that it is indeed possible in this framework 
to separate the chiral modes. To be able to do so analytically, we now let 
L s — > oo, and we switch to Kaplan's version, with an ordinary Wilson action 
in the s direction plus a mass-like term m(s) = m ■ sign(s), sign(s) = 1 for 
s > 0, —1 for s < 0, and at s = 0. Hence the free fermion action reads 



1 



_l — Y^^ s s , +l _|_ S s+l s , - 25 S}S >)]6 r fi 



^x+r,s' , (24) 



where p^, A refer to the HF at some mass M. We now follow the procedure 
of Ref. [34] to look for solutions of the type ^_r,l = ^ px ^sUR,L > where 



3 Generally, the changes M — > — M and a — > —a imply Pn — » Pn, A 
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15Ur = ur, 75UL = —ul. We require the inverse propagator G 1 (p,s) to 
reduce to a chirally invariant G~ 1 (p), which amounts to the condition 



^(S s ^ + i-5 s+ i^)^+m(s)-^^-(5 S)S ' + i+5 s+ i yS '-25 SyS >)+X(p) 



<S> s ur,l = 0. 
(25) 

Following Kaplan [[H]] we make the ansatz = z$ s for s, s + 1 ^ 0. The 
solutions for z are 



*7± JU 2 -r^ s + l 

^?,L = , U = m(s) + A(p + r w>s . (26) 

rw,s T 1 

The upper (lower) sign in the denominator refers to zr (zl), and in both 
cases the two signs in the numerator are possible. As an example, in the 
limit ryy tS — > 1 there is only one finite solution: zr = 1/U. Normalizability 
requires 

\z R (s>0)\ = |- 1 | < 1 -> \{p)>-m or A(p)<-(2 + m) 



1 + m + A(p) 
1 

1 — m + X(p) 



\ ZR (s < 0)| = |- 1 | > 1 -> (m - 2) < A(p) < m . (27) 



The allowed region for A(p) is shown as shaded areas in Fig. || on top. 
Below we show the interval of values that X(p) actually takes, depending on 
M (for a = aj\j)- The minimum is at X(p = 0) = M 2 /(e M — 1), and we 
see that there are many possible combinations of m and M which yield a 
single right-handed fermion. In some cases the solution collapses at some 
"critical momenta", in others it extends all over the Brillouin zone. It is 
particularly favorable to choose a rather large parameter M, so that 4d 
locality is excellent and - as a related property - X(p) is confined to a 
narrow interval for all momenta. Then, for instance m ~ 1 ... 2 guarantees 
the solution to exist at any momentum. 

In the limit ryy,s — > — 1 the situation is analogous, but there one deals 
with a left-handed fermion. 



IV. Conclusions 



We have revisited the fixed point action approach to the formulation 
of lattice chiral fermions, emphasizing in particular the blocking from the 
continuum and the role of the Ginsparg- Wilson relation. We also consider 
the artifacts in the remnant chiral symmetry in short ranged truncated fixed 
point actions. Based on a formula used in the overlap formalism, we found 

14 Here we mean full chiral invariance, not just the GWR. Hence G~ 1 (p) is not a trun- 
cated perfect Dirac operator, but just its vector part «p M (p)7 M , which - in its isolated form 
- would imply fermionic doubling. 
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a method how a lattice Dirac operator can be "GWR corrected" . This leads 
to a large class of new solutions, and it allowed us to further optimized 
their locality, which pays off in reduced spectral artifacts after truncation. 
Along these lines, we also arrived at a formulation of improved domain wall 
fermions. We discussed practical issues in view of their application in simu- 
lations of vector-like theories. The improvement could, for instance, reduce 
the extent L s in the fifth direction, which is needed for a number of favorable 
properties. Finally, we showed that the mechanism to separate the chiral 
fermions still works in our improved formulation. 

Acknowledgment I am indebted to S. Chandrasekharan, F. Niedermayer, 
K. Orginos, M. Peardon and U.-J. Wiese for interesting comments. 

Appendix 

We know that FPAs obey the GWR, which implies plenty of attractive 
properties, but we also know that such actions are extremely hard to im- 
plement. What we need is a good but very simple gauging of our HF, and 
here we want to comment briefly on the possibility to use the GWR for the 
construction of such a gauging. 

We denote the lattice Dirac operator by D and refer to a = 1/2, so that 
the GWR can be written as 

D{r) +D(r) =^D{x)D{r-x) , D = 75D75 . (28) 

X 

Let D{r) represent the HF in the free case, which obeys the GWR to a good 
accuracy (see Table |l|). The task is now to find a gauging which preserves 
this property in the interacting case. This is rather involved, because now 
the GWR is supposed to hold approximately for each lattice path. (The 
situation is even worse in cases where a 7^ const.). 

As a simple example, we consider the gauging by the average of the 
shortest lattice paths only, and we discuss the case d = 2. Regarding con- 
dition (^), the case r = (00) remains exact, and for r = (11), (21) and (22) 
the combination of couplings, which nearly matches in the free case, acts on 
all paths involved, hence the GWR is still approximated well. For r = (20) 
that combination of couplings splits into coefficients for the paths of length 
2 and 4, but all those coefficients are small, so the situation is not much 
worse than in the free case. The main problem is the case r = (10), which 
amounts to 

- 0.490 [ link } = -0.608 [ link } + 0.061 [ sum over staples } . (29) 

Here the matching path by path does not work well, and this provides some 
insight into the limitation of the minimal gauging by shortest lattice paths. 
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One should now correct for that by inserting fat links, a clover term etc. 
and tuning its coefficients. This work - and the extension to d = 4 - is in 
preparation. 

At this point, we just add a remark on how to simplify this task. If we 
generally define X = 1 — D, X = 75X75, then the GWR for a = 1/2 reads 

11 = 1. (30) 

If D is linear in the 7^'s, D = p^j^ + A, then the GWR is identical to the 
requirement that X be unitary. We insert again the HF for D and introduce 
a new hypercubic variable A = 1 — A, so the GWR takes the form 

Y^~P^P^ r ~ x ) + AO)A(r - x)) = 5 rfi . (31) 

X 

Again this holds exactly at r = for any reasonable gauging, but r ^ is 
complicated. Of course, we cannot avoid the inconvenient quadratic term 
by these re-definitions, but at least the rest has been optimally simplified. 
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Figure 1: Spectrum of the "hypercube fermion" (on top) and its GWR 
corrected and truncated modification (below) in d = 2. In both cases, the 
spectra keep close to the unit circle, hence the artifacts due to truncation 
are small. In the most contaminated region around e ~ 2 we observe some 
progress thanks to GWR correction. 



19 



CO 

e 




05 1 1.5 2 2.5 3 



Re 8 



OJ 

I 

O.E 







-0.5 



-1 




0.5 1 15 2 2.5 



Re 8 

Figure 2: Spectrum of the "hypercube fermion" (on top) and its GWR 
corrected and truncated modification (below) in d = 4. Again, the spectra 
keep close to the unit circle. From the shape the progress due to GWR 
correction is not easily visible here, but if we measure the mean deviation 
from the unit circle, which takes into account the eigenvalue density, then 
we observe an improvement by 31 percent. 
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Figure 3: On top, the shaded regions are the allowed areas for X(p) providing 
a right-handed Kaplan-type fermion for rw,s = 1 an d m(s) = m • sign(s). 
Below, we show the interval where X(p) takes its values for a given parameter 
M and a = olm- (For a = —a-M the shaded area is mirrored at the origin.) 
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